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ABSTRACT 

Aims. A study of the time sequence of optical colours of a large sample of nearby Type la supernovae (SNe la) is presented. We study 
the dependence of the colour time evolution with respect to the lightcurve shape, parametrized by the stretch factor. 
Methods. We fit the spectral template that minimizes the colour dispersion in SNe la, as measured through UBVRI photometry of 
near-by supernovae. 

Results. A clear colour dependence upon lightcurve shape is found, with the narrower lightcurves being redder up to about one month 
past lightcurve maximum. We also derive an average reddening law, after correcting for lightcurve shape differences in intrinsic colour, 
that is well described by a Cardelli, Clayton & Mathis law with R v = 1.75 ± 0.27 for 80 Type la supernovae with E{B — V) < 0.7 
mag. A subset sample including 69 SNe with modest reddening, E(B - V) < 0.25 mag, yields a significantly smaller value, R v ~ 1, 
suggesting that the observed reddening of Type la supernovae may have a more complex origin, perhaps involving other processes 
beside extinction by interstellar dust in the host galaxy. 

Key words, supernovae: general - Stars: statistics 



1. Introduction 

The use of type la supernovae (SNe la) as standardized can- 
dles has been successfully exploited in recent years as a precise 
tool for cosmology to discover the accelerati on of the Universe 
through some yet unde termined dark energy dRiess et al.lll998t 
iPerlmutter et al.l U999I) . More recently we have seen the start 
of the second generation of supernova surveys, with dedicated 
supernova searches (e.g SNLSQ SDSS0, ESSENCE^) and mas- 
sive spectroscopic follow up. The results from these surveys are 
showing that increasing the number of supernovae populating 
the Hubble diagram at high redshift may not be sufficient to dis- 
criminate among alternative dark energy models. Systematic un- 
certainties appe ar to be the limiting factors for the ongoin g su- 
pernova efforts (lAstier et al.ll2006t IWood- Vasev et aill2007h . 

Dimming by dust along the line-of-sight, predominantly in 
the host galaxy of the supernova explosio n, is one the mai n 
sources of systematic uncertainty, see e.g. iKnop et all (|2003). 
For Type la supernovae, two additional causes of reddening are 
sometimes suggested: extinction by circumstellar dust and the 
possibility of an intrinsic colour-brightness relation. The last two 
scenarios further complicate the prospects to disentangle the ef- 
fects and their possible evolution with cosmic time. 

Our ability to correct for these effects relies on the knowl- 
edge of the intrinsic colours of SNe la, as well as on the nature 
of the reddening. The latter is often assumed to be similar to 
what is found for extinction by dust in the Milky Way. As we 
will see in section [6] this is not necessarily the case, and needs 
further investigation. 
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The effect of cosmological redshift on the measured broad- 
band filters is accounted in a s traight-forward manner through 
^-corrections (Ki m et al]|l996l) . However, the SN spectrum is 
usually not measured at all epochs, since typically only few 
spectra are taken to allow type identification. The standard tech- 
nique consists in using spectroscopic templates, built by aver- 
aging spectra of well observed (mostly nearby) SNe la. Thus, 
the uncertainty in jT-correctio ns depends pr i marily on the spec- 
troscopic diversity of SNe la. iNugent et al.l d2002l) showed that 
changes in individual spectral features do no significantly af- 
fect ^-corrections, as long as there is a good match between 
the observed and the rest frame bandpass filters. These are in- 
stead mainly sensitive to the supernova colours. Once again, the 
knowledge of intrinsic SN colour can help to solve one of the 
limiting factors in supernova cosmology. 

This issue becomes even more critical when dealing with the 
observations of very high-z SNe, for which the rest frame ex- 
tends, partly or entirely, into the UV part of the spectrum. Our 
limited knowledge of the supernova properties in the U-band, 
due to poor telescope and CCD sensitivity at these wavelengths, 
was identified as the main source of uncertainty i n the determi- 
nation of cosmological parameters in lKnop et al.l d2003l) . 

Some efforts to observe SNe in the U-band have been made 
and collected a number of nearby SNe da ta dJha et al.l 2006; 
Pastorello et al. 2007; StanishevetaD|2007). At the same time 
there have been attempts of using spectra of high redshift SN la 
from the SNLS sample, redshifted to optical wavelengths with 
good sensitivity, under the as sumption that there is no evolution- 
ary tr end in the SN spectrum dGuv et al.l2005l 120071; fHsiao et al.l 
2007). This assumption may not hold up to arbitrary redshifts. 
As the average metallicity of the universe increases with cos- 
mic time, it is not unreasonable to expect that high redshift 
SNe la stem from environments with lower average metallicity. 
The effect on the spectral energy distributi on of a lower m etal- 
licity progenitor has been modeled bv lHoeflich et al.l (1 19981) and 
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iLentz et all (|2000) who found that such SNe la, especially at 
early epochs, are expected to show enhanced flux in the UV re- 
gion of the spectrum, weaker absorption features in the optical 
and a shift in the minima of optical absorptio n features to redder 
wavelengths (see, however. lEllis et all (120071) ). 

In this paper we present a statistical study of U - B, B - V, 
V - R and R — I colour evolution for 80 nearby supe rnovae. This 
work i s a continuation of the analysis presented in Nobil i et ail 
(120031) . with an increased sample and with the addition of U- 
band data. This allowed us to study the dependence of colours on 
the stretch factor. That is, in turn, applied to correct the spectral 
templates used for computing the ZT-corrections. The knowledge 
of intrinsic colour dispersion is further used to study the average 
properties of the extragalactic reddening law of SNe la. 

In section [2] the data set used in this paper is presented. 
In section [3] we investigate the correlation of colours with 
lightcurve shape, and estimate colour curve models for typical 
SNe la. The study of intrinsic dispersion in colours is reported 
section |4] report. In section [5] we discuss the use of our colour- 
stretch relation to modify the spectral template for computing 
Zf-correction, and compare our results with previously published 
estimates. Finally in section[6] we discuss the implication of our 
analysis for computing the effective average extinction law of 
SNe la. 



2. The data set 

We present the analysis of the colour properties of a collection 
of 80 well observed nearby SNe la availa ble from the litera- 
ture. This includes th e Calan/Tololo da ta set dHamuv et al.l l996) 
and th e CfA data set (iRiess et al.ll 1998b . an alyzed in l Nobili et al.1 
(120031) . and the more recent data set by iJha et al.l d2006l) . to- 
gether with some single very well observed nearby supernovae 
(SN 200 lei (iKrisciunas et aTl l2003l). S N 2000ca and SN 200 lb a 
jKrisciunas et al.ll2004a[): SN 200 lcn dKrisciunas et al.ll20 04b): 
SN 2001V dVinko et alfeO Ol; SN 2002er, dPignata et al.ll2004l): 
SN 2003du (IStanishev et al.ll2007l) : SN 2005cf dPastorello et al.1 
120071) ). Sub-luminous Type la SNe, 1991bg-like, were not in- 
cluded in the sample nor supernovae with poor Z?-band sampling 
around lightcurve maximum. Out of these 80 SNe, a subset with 
moderate colour excess, E(B -V)< 0.25 mafl was selected to 
study the intrinsic colours of Type la supernovae. The complete 
list of the SNe used is reported in Table Q] together with the ob- 
served filter data available, redshift and fitted g-band s tretch fac- 
tor (s) dPerlmutter et alJll997t iGoldhaber et al.ll200ll) . Figure Q] 
shows the distribution of stretch factors and redshifts for the 
whole sample. 

The d ata have been A"- corrected to the rest frame bandpass, 
following Kim et al.l d 19961) . assuming Bessell filter transmission 
curves dBesselllll990l) . The spectra l template u s ed for computing 
the ^-corrections was obtained by Hsi ao et all d2007l) by averag- 
ing about 600 spectra of SNe la at different epochs. The sample 
used for this analysis includes nearby supernovae with redshifts 
up to -0.1, with most of the SNe at z < 0.06 (see Fig. [T}. Our 
ability to correct for extinction, i.e. to compute the colour excess 
E(B - V), depends on our knowledge of the intrinsic colours of 
SNe la. Thus, an iterative procedure has been adopted until the 
procedure converged, as explained in the next section. 
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Fig. 1 Distribution of the B-band stretch factor, s, (top panel), 
distribution of the redshift (bottom panel) for the low extincted 
sample of 69 SNe. 

3. The colour-stretch relation and colour time 
evolution 

The correlation between colours and lightcurve shap e parame- 
ter is a known property of SNe la. Phillips et al. d 19991) showed 
the dependence of B — V and V — I at lightcurve maximum on 
the shape parameter. We use the larger data set available to ex- 
tend the analysis to all epochs up to 60 days after lightcurve 
maximum. We also expand the wavelen gth range to also include 
the t/-band. As in iNobili et all d2003l) . the SN colours are de- 
rived from the data points, without any interpolation. For this 
to be possible, the SN should be observed in at least two dif- 
ferent bands at the same epoch, which is often the case. Thus, 
the Z?-band lightcurve fit is performed exclusively with the aim 
of determining the time of maximum and the stretch factor. We 
consider the rest frame time since B-band maximum, t, and the 
stretch factor, s, as two independent variables. We also introduce 
the colour excess E(B - V) for each SN as independent variable. 
Each colour U-B, B-V, V-R and R - I is fitted with a 
parametrised function: 



X - Y = b XY (t) + a XY (t) ■ (s - 1) + c XY ■ E(B - V) 



(1) 



4 In section 6 we will relax this constrain, adding the remaining 1 1 
SNe, and study the impact on the derived reddening law 



where X and Y are arbitrary filters, a(t) and b(t) are determined 
as least squares cubic spline fit of 6 knots each, distributed in the 
time range between -10 and +60 days with respect to the time 
of maximum light. The knot positions are chosen to minimize 
the number of knots while yielding a good fit to the data. Note 
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that, in this parametrisation, the b functions correspond to the 
colour curves for a stretch factor 5=1 SN, while the a functions 
give the dependence of the colours on stretch. Different depen- 
dences on s t retch h ave been investigated, e.g. (1/s 3 - 1) used by 
Knop et all (120031) . However, the simplest choice of Eq. Q] de- 
scribes the current data best. The cxy parameters are fitted for 
each of the colours without assuming any specific wavelength 
dependence of the dimming of Type la SNe. 
An interactive procedure was applied as follows: 

1 . Fit the B-band lightc urve time of maximum and the stretch 
parameter, following Goldhab er et al.l (1200 ll) 

2. Compute the colour U — B, B— V, V — R and R — I for each 
SN at each observed epoch 

3. Correct the observed colour s for Milky- Way ext inction, as- 
suming the extinction law (ICardelli et al.| [l989)(CCM). as 
modified by |0 'Donnelll d 19941) . The corrections are com- 
puted by reddening the SNe la spectral template at different 
epochs, and integrating in all band-passes. 

4. Compute /f-corrections for each SN using the spectral tem- 
plate corrected by the colour-stretch relation determined in 
the previous iteration. In the first iteration, Zf-corrections 
where computed without using any colour-stretch relation, 
i.e. as if all SNe had s = 1. 

5. Compute the colour excess E(B - V) as weighted average 
excess for all epochs from the B - V colour curve for each 
supernova determined in the previous iteration. The inverse 
square of the uncertainty on the individual colour measure- 
ments, 1/cr?, is used as weights in the average, and the uncer- 
tainty on E(B- V) is determined as 1 / Jj^l/o?. For the first 

iteration, an initial colour excess E(B - V) for each SN was 
estima ted following the method described by Phillip set alJ 
(119991) . 

6. Fit the function in Eq.[T]to the data separately for each colour. 

7. Modify the spectral template for the average colour-stretch 
relation. 

The spectral template colours are modified by using a cubic 
spline interpolation of the ratio between the synthetic photom- 
etry and the new photometry at the effective wavelengths, for 
the U, B, V, R and /-band. These are determined from the fitted 
colour curves as follows: 

U = (U-B) + B' 
B - B' 

V = B'-(B-V) (2) 
R = V - (V - R) 
I = R-(R-I) 

where B' indicates the B-band lightcurve by iGoldhabe r et al.l 
(1200 ll) . and the colours between brackets are the one fitted in 
the previous iteration. 

The steps 1-7 were repeated until the difference in the colour- 
curves between consecutive iterations changed by less than 
0.5%. Only three iterations were needed to reach such precision. 

Tables [5]- [8] report the values for a and b for each colour as 
a function of time. Figure [2] shows the time evolution of U - B, 
B - V, V - R and R — I for different values of the stretch factor, 
s = 0.8,0.9, 1.0, 1.1, 1.2. Figure [3] gives a comparison of the a 
parameter fitted for each of the colours. For U - B the depen- 
dence on stretch is substantial already before maximum light, 
and decreases with time (see discussion in Section[8j. In all other 
cases, the dependence on stretch peaks between 15 and 20 days 
after fi-band maximum. This coincides with the epoch at which 



the photosphere recedes into the core of the supernova, and ab- 
sorpti on features from iron gr oup elements appear in the spec- 
trum. iKasen & Wooslevl (120071) show that the colour evolution 
depends on the ionization evolution of the iron group elements, 
which is faster for dimmer supernovae. Their model dependence 
is in good agreement with the results of this work. 

Note that the data points were considered to be uncorrelated 
at different epochs when performing the fit. Thus, even if some 
SN lightcurve are better sampled than others, all points have the 
same weight in the fit. As this is only a crude approximation, we 
take this into account in the determination of the uncertainties. 
For this reason, the uncertainties on the colour curves reported 
in Tables |5]-|8] are determined by the use of Monte-Carlo simu- 
lations, as explained in Section [7] 

4. Intrinsic colour dispersion 

The found colour-stretch relation contributes to decreasing the 
dispersion along the colour curve for all colours, once it is ap- 
plied to the data. However, the remaining dispersion is still too 
large to be explained solely by statistical fluctuations based on 
the measurement errors. This forces us to conclude that there 
is remaining intrinsic di spersion in the SN colours. We follow 
the same procedure as in lNobili et al.ld2003l) in order to estimate 
the intrinsic dispersion for each colour as a function of time. 
The data, corrected for the colour-stretch relation found in the 
previous section, are divided in time bins. For each time bin, we 
compute the weighted average of the residual to the colour curve 
(for a s — 1 supernova), and the sample weighted standard de- 
viation, as the squa re root of the weighted second moment (see 
iNobili et all d2003l) ). Tabled gives the results of this analysis. 

The weighted standard deviation can be taken as an estimate 
of the intrinsic disp ersion in each time bin, A, however, as al- 
ready pointed out in lNobili et al.1 (120031) . this is an overestimate 
since part of the scatter is due to measurement uncertainty only. 
In order to disentangle the two contributions, we run a Monte 
Carlo simulation to generate synthetic data sets, with dispersion 
given by the quoted measurement uncertainties only. Thus, we 
compute the weighed standard deviation on the simulated data 
sets, 6. Finally we run a hypothesis test, to compare the disper- 
sion measured on the real data to those measured on the synthetic 
data sets for each colours. 

We set a level of significance a = 0.01 for rejecting the null 
hypothesis, i.e. the probability that A + 6 under the assump- 
tion that the null hypothesis is co rrect is assumed equal to 1% 
(INobili et al.l 12003b ICowanlll998l) . For all cases for which the 
null hypothesis is rejected, we compute the quantity: 

A con " = VA 2 - 6 2 (3) 

that we take as an estimate of the intrinsic dispersion, and 
a lower limit on this value is set at a 99% confidence limit. 
We found no cases compatible with no intrinsic dispersion, (see 
Tabled). 

5. Supernova spectral template 

The spectral template built by iHsiao et all (120071) . ("hsiao" in 
our plots), was used as the starting assumption in the iterative 
method for this analysis. The template was built using ~600 
spectra for ~ 100 SNe la where the observed spectra were colour- 
corrected to match a single SED, before averaging them into 
one uniform spectral template. The colours used were those in 
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Fig. 2 From top to bottom, U — B, B — V,V — R and R - I time evolution for different stretch values: s = 0.8, 0.9, 1.0, 1.1, 1.2. The 
straight line at late time in the B - V evolution is the Lira line. The abscissa is the restframe epoch since the B-band maximum light. 



Kno p et ail d2003l) . corresponding to a "normal" 5 = 1 SNe la. 
The hsiao spectral template is meant to be used for computing K- 
corrections after correcting it for the observed colours. This has 
the advantage of being independent of an assumed stretch-colour 
relation, but the resulting uncertainties on the A"-corrections, de- 
pend on the availability of observed colours, and can b e as large 
as 0.12 mag, as showed in Fig. 9 of iHsiao et alj d2007l) . By us- 
ing the colour-stretch relation derived with our technique, the 
spectral template colours can be adjusted for each supernova 
lightcurve shape, before computing -corrections. In our case, 
the uncertainties on A"-corrections are dominated by the intrin- 
sic dispersion in the supernova colours. Whether these corre- 
spond to an intrinsic variability of the spectral features is be- 
yond the scope of this work. We note, however, that there are 
several indications that at least some of the characteristic SN la 
spectral features correlate well with lightcurve shape parame ters 
dNugent et alJll995HBenetti et al.ll2005l:lGaravini et alll2007h . 

Another interesting approach is the one in salt2 dGuv et alj 

|2007|) . Instead of building spectral templates to be used for com- 
puting A"-corrections, an empirical model that describes the time 
variation of the spectral energy distributions and its dependence 
on a lightcurve shape parameter (corresponding to stretch) is 
used. In this case, the SED is a function of epoch, wavelength 
and lightcurve shape parameter. This ambitious model is trained 
on about the same data set as the hsiao templates. However, due 
to the larger amount of parameters needed to describe the model, 
the uncertainties are quite large, especially in the U and / part of 
the spectrum, where fewer measurements are available. 



A comparison between the spectral templates for a s = 1 su- 
pernova at some epochs from maximum light (day=-10,0,10,20) 
is shown in Fig. [5] Differences in our colour corrected spectral 
template and the original hsiao template are small, up to 2%, 
and even smaller around maximum, about 0.5%. This differ- 
ences originate from different data sets used to deter mine intrin- 
sic co lours. The HsiAO-template uses the colours in Kno p et ail 
(I2003h . based on a smaller sample than the one used here. We 
note in particular, the U - B colour that was derived based on 
measurements of only 5 SNe la at maximum. When comparing 
our spectral template with the salt2 template we find large dif- 
ferences already for s = 1 . We note, however, that the differences 
are mainly in the U part of the spectra where large uncertainties 
of the salt2 templates are reported by the authors. 

To further compare the spectral templates, we have com- 
puted A"-corrections using salt2 templates, hsiao templates and 
our templates for different values of the redshift. The results of 
the comparison for 5=1 are reported in Table [3] The larger 
discrepancy is found for z = 0.6, where the mismatch of the ob- 
served and restframe filter is largest. Th is is in agreement with 
the results found by IHsiao et ail d2007l) . suggesting an uncer- 
tainty of ~ 0.04 in A"-corrections for this redshift. 

A comparison between the dependence of A"-corrections on 
stretch, is only possible between salt2 and our templates. Fig. [6] 
shows the A"-corrections to restframe B-band for two different 
redshifts z = 0.3 and z = 0.6 computed using the two sets of 
templates. The r.m.s. on the difference is reported in Table[4] 
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Fig. 3 From top to bottom, the a parameter time evolution for U — B, B— V, V — R and R - 1 (solid line), and its uncertainty (dotted 
line). The minimum in the curves from top to bottom occurs at day 15, 17, 20 and 20. 



All the differences found demostrate that by using the un- 
corrected spectral template, it is possible to introduce systematic 
uncertainties as large as 0.04 mag, limiting our ability to use 
SNe la for precision cosmology. 



6. The reddening law from SNe la 

Corrections due to the wavelength dependence of the dimming 
of the supernova light is one of the largest systematic uncer- 
tainties in supernova cosmology. While the empirical evidence 
of the dependence of peak magnitudes on SN colours is very 
clear, there is no consensus on how to disentangle the contri- 
bution from reddening due to dust in the host galaxy and in- 
trinsic SN colour variations. Moreover, most published results 
assume the dust in the host galaxy to have similar properties 
as (the average) Milky Way dust, implying a total to selective 
extinction parameter Rv = 3.1. Studies of dust properties in 
small sets of distant galaxies indicate mean values compatible 
with this value, even though a significant range of values may be 
prese nt (lOstman et al J 120071: iGoudfrooij et all 1 1994 iPatil et al.1 
120071) . Ideally, cosmological distance estimates using SNe la 
should include extinction corrections for each particular line of 
le cific dust propertie s 
dNordin et alj|2008b . 
between lines of sight are much less problematic for cosmologi- 
cal applications than a potential bias in the adopted mean value. 

Intriguingly, there are a large number of examples of mea- 
surements of the mean Ry from SNe la, based on different meth- 
ods, that disagree with the universally assumed value: a sig- 



sight, including the spec ific dust propertie s of each host galaxy. 
However, as shown in (Nordin et al. 2008), the variations of Ry 



nificantely smaller valu e of the mean Ry is f ound. Some ex- 
ampl es are Ry ~ 1.8 dKrisciunas et al. l l200Cl). Rv = 2.55 ± 
0.30 dRiess et al.1 1 19961). R v = 2.6 + 0.4 dPhillips et al l [l999). 
R v = 1.09 (iTripnl 119981) . Ry = 2.5 (Altav iHaet allEiol 
B = 1 77 ± 0.16 corres ponding to R v = 0.77 + 0.16 dAstier et all 
120061 iGuv etani2007l) . 

Using the fitted c parameters in Eq. Q]the reddening law that 
best describes the data can be determined. In fitting the c param- 
eter for each colour, cub, cbv, cvr, Cri, we minimize the part of 
the colour dispersion that depends on reddening, without assum- 
ing any value for Ry. The four c parameters, one for each colour 
fitted, can then be used to estimate the parameter Ry : 



cub = Ru ~ Rb 
cbv — Rb — Rv 
cvr — Ry —Rr 
cri — Rr- Ri 



(4) 



where Rx = Ax/E(B - V), and Ax is the dust absorption in a 
given band. We note that cbv = 1 by construction (see EqQ~|i, 
so we have 3 independent equations. Fig. [8] shows - As as 
a function of wavelength for E(B - V) = 0.1. The data points 
are determined using the c parameters. The dotted line is the 
Cardelli, Claython & Mathis (CCM) law assuming R v = 3.1 
and the light grey solid line is the result of fitting the CCM law 
on our data points. The best fit is obtained for Ry = 1.01 ± 0.25. 
The dashed line is the result of a similar procedure obtained by 
dGuv et al.ll2007l) while developing the salt2 templates. A typi- 
ca l value of E(B — V) = 0.1 was used in analogy with the work 
of iGuv et al.l (12007). We note that our measurements agree well 
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Fig. 4 From top to bottom, left to right, U — B, B— V, V — R and R - I time evolution compared with the data corrected for Milky 
Way extinction as well as for the colour-stretch relation and host galaxy extinction following Eq. [1] 



with the results from salt2 in the R and /-bands, but disagree 
in the <7-band. We note also, that the salt2 curve is not easily 
modeled by simply assuming a different value for Rv in CCM. 

For testing the effect of selecting a low extincted sample, 
we have introduced in the analysis 1 1 high extincted SNe. We 
assumed the same colour curves obtained on the low extincted 
sample, and we investigated the reddening law on the larger sam- 
ple. The c parameters obtained in this way are slightly different, 
and lead to a larger value of Ry when fitted by the CCM law, 
equal to R v = 1.75 + 0.27. Given that the two samples are not 
independent from each other, since 69 SNe are present in both 
samples, the different values found are strikingly incompatible. 



7. Monte Carlo simulations 

Due to the potential correlation of the fitted parameters a, b and 
c, and the possible correlation in the intrinsic dispersion, we have 



chosen to use Monte Carlo simulations to estimate the uncertain- 
ties in the above parameters, as well as on Rv- 

The data suggests that the intrinsic colour dispersion is cor- 
related between points at different epochs, e.g. if a SN is bluer 
than average at a given epoch, it is likely to remain bluer0 

However, the correlation coefficient as a function of epoch is 
quite difficult to quantify, given the small sample and the sizable 
measurement uncertainties. Moreover, the inhomogeneity of the 
data sample could partially introduce such correlations, espe- 
cially in the f/-band, for which a false correlation could originate 
by problems in the S -corrections between filter systems (see also 
the discussion in the next section). For these reasons, when run- 



5 We note that, the colour excess E(B - V) determined from the B — V 
colour, is used for the fit of Eq. [T]to all the other colours. Thus, a SN 
can still remain redder or bluer than average even after the reddening is 
taken into account. The reddening law that minimizes the scatter around 
all the colours is the result of the method used. 




Fig. 5 Comparison of the spectral templ ates derived in this work for s — 1 (black solid line), the spectral template by iGuy et al.l 
(2007) (blue dashed line) and the one by Hsiao et al. (2007) (red dashed-dotted line) for epochs=-10,0,10,20referred to the time of 
maximum. 



ning the Monte Carlo simulations, we considered the uncertainty 
in the data for the two extreme cases of either full correlation or 
no correlation at all in the intrinsic dispersion (see Section 0J. 
We found the case with no correlation gives the larger uncertain- 
ties on the fitted parameters, and we have conservatively chosen 
these as our estimate of the uncertainties. 

We generated 100 synthetic data samples, with as many 
lightcurves as in the real set, according to the following prescrip- 
tion: 

- We generate data points normally distributed around the light 
curves, with a dispersion, cr, given by the measurement un- 



certainties, assumed to be uncorrelated, plus an intrinsic dis- 
persion fully correlated in time, i.e. all points of the same 
SN for a given band where given the same intrinsic disper- 
sion at all epochs (case A) and fully uncorrelated (case B). 
The intrinsic dispersion is considered Gaussian, with a stan- 
dard deviation of 0.05/ V20. The data points are generated 
at the same epochs as the real data. 



6 The value of 0.05 has been chosen a posteriori as average intrinsic 
dispersion in all colours at all epochs, given the results shown in Table[2] 
Using a smaller value, such as 0.03, closer to the lower limit values in 
Table 12 we find smaller uncertainties on the parameters. In addition, 
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Table 2 Results of the analysis of all SNe. First column: central 
value in days for each time bin; Nk is the number of points for 
each bin; Axy is the intrinsic dispersion computed as weighted 
standard deviation; At™ T is the corrected intrinsic dispersion, 
computed as in Eq. [3] and in the last column is the estimated 
lower limit at 99% C.L. (see text); 



day 


N k 


Aub 


ACOIT 


L.L. 


-10 


22 


0.13 ± 


0.03 


0.13 


0.09 


U 


Q1 

a J 


0.10 ± 


0.03 


(1 in 

U. 1U 


a aa 
U.Uo 


10 


67 


0.06 ± 


0.05 


0.06 


0.03 


20 


44 


0.10 ± 


0.02 


0.09 


0.06 


30 


29 


0.06 ± 


0.07 


0.06 


0.04 


40 


20 


0.07 ± 


0.02 


0.06 


0.05 


day 


N k 




A can - 


L.L. 


-5 


77 


0.05 ± 


0.01 


0.05 


0.03 





180 


0.03 + 


0.01 


0.03 


0.02 


1 A 

IU 


TOO 

22s 


0.08 ± 


0.01 


a no 
U.Uo 


0.04 


1 9. 
1 ft 




0.05 ± 


0.01 


n ac 


U.UZ 


30 


137 


0.10 ± 


0.01 


0.09 


0.05 


40 


74 


0.06 ± 


0.01 


0.06 


0.03 


50 


66 


0.12 + 


0.04 


0.12 


0.06 


day 


N k 




ACOIT 


L.L. 


-5 


64 


0.04 ± 


0.01 


0.04 


0.02 





120 


0.04 ± 


0.01 


0.04 


0.02 


10 


131 


0.03 + 


0.01 


0.03 


0.02 


15 


65 


0.08 ± 


0.01 


0.08 


0.04 


20 


82 


0.03 ± 


0.01 


0.03 


0.01 


30 


81 


0.03 ± 


0.05 


0.03 


0.02 


40 


58 


0.05 ± 


0.01 


0.05 


0.03 


50 


47 


0.03 ± 


0.01 


0.03 


0.02 


day 


N k 




ACOIT 


L.L. 


-10 


44 


0.04 ± 


0.01 


0.04 


0.03 





147 


0.06 ± 


0.01 


0.06 


0.03 


10 


141 


0.10 ± 


0.02 


0.10 


0.05 


20 


77 


0.13 ± 


0.02 


0.13 


0.07 


25 


97 


0.05 ± 


0.02 


0.05 


0.03 


40 


62 


0.14 ± 


0.03 


0.14 


0.08 


50 


50 


0.14 ± 


0.04 


0.14 


0.08 



- We compute the colours and use the measured E(B - V) to 
add reddening in each colour, given an assumed value of Ry. 

- We fit the colour curves for the simulated samples, and re- 
trieve the value of Ry using the fitted c parameters. 

The dispersion in the fitted parameters in all the simulated 
data sets is taken as uncertainty in the parameters, a, b and c 
fitted on the real data, as reported in Tables [5]- [8] 

We run this test for two values: Ry me) - 1 and Ry me) = 1.75. 
In both cases we were able to retrieved the true value, i.e. we 
found no bias in our method. The dispersion measured as r.m.s. 
on the distribution of the Ry can be taken as an estimated upper 
limit to the uncertainties on the determination of Ry. We found 
<x = 0.13 (case A) and <x = 0.27 (case B). 

8. Discussion 

An extended analysis on supernova colour-curves and their de- 
pendence on the light curve shape parameter was presented in 
previous sections. An empirical model was defined, described 
by Eq. Q] and it was succefully fitted to the data. The robustness 

we assumed the intrinsic dispersion in colours to be equally distributed 
in the two bands 



Table 3 Comparison at different redshifts between the K- 
corrections computed using the various spectral templates for 
s — 1. The corrections are computed from the band listed in 
the first column to restframe B-band. The mean and the r.m.s. 
are computed over time on the residuals to the ^-corrections 
computed using our templates. The largest dispersion is noted 
at z = 0.6 where the mismatch betw een observed and restframe 
filters is maximum, as also noted bv lHsiao et al.l (120071) 



Band 


z 


mean 


r.m.s. 


HSIAO 


B 


0.1 


0.008 


0.020 


V 


0.3 


0.007 


0.015 


R 


0.6 


0.012 


0.031 


I 


0.8 


-0.011 


0.028 


I 


0.9 


0.005 


0.010 


salt2 


B 


0.1 


0.023 


0.017 


V 


0.3 


0.021 


0.015 


R 


0.6 


0.041 


0.040 


I 


0.8 


-0.040 


0.026 


I 


0.9 


0.016 


0.013 



Table 4 Comparison between the dependence of Zf-corrections 
on stretch for two different redshifts, for salt2 and our templates. 
The mean and the r.m.s. are computed over time on the residuals 
to our templates. Once again the largest dispersion is noted at 
Z = 0.6 where the mismatch between observed and restframe 
filters is maximum, and for more extreme stretch values. 



z 


s 


mean 


r.m.s. 


0.3 


0.8 


0.014 


0.045 


0.3 


0.9 


0.019 


0.028 


0.3 


1.0 


0.022 


0.015 


0.3 


1.1 


0.023 


0.013 


0.6 


0.8 


0.061 


0.123 


0.6 


0.9 


0.049 


0.078 


0.6 


1.0 


0.042 


0.041 


0.6 


1.1 


0.036 


0.033 



of the analysis, was thoroughly tested. The residuals from the 
colour curves were analyzed to search for further dependence on 
the stretch factor, SN redshift, or colour excess E(B - V). No sig- 
nificant correlation was found. Furthermore, the colour curves 
we derive describe well the subsample of very low extinction 
objects, E(B - V) < 0.05, suggesting that the analysis is robust 
and not particularly affected by moderately reddened SNe. 

In general, the dispersion we measured in all colours at all 
epochs is non-negligible, indicating that colours of Type la su- 
pernovae have more scatter than can be accounted for by the 
measurement uncertainties. These findings do not support ear- 
lier results based on a smaller sample where a neglig ible intrin- 
sic scatter at late epochs was found in some colours (Nobili et al. 
120031) . The same conclusion was reached by iJha et al.l (2007) 
when studying an overlapping sample of nearby SNe la. They 
found an intrinsic dispersion of cr B _ v = 0.062 at epoch +35, 
in agreement with the result reported here. The intrinsic disper- 
sion for various epochs and colours is reported in Table [2] and 
could potentially be used to constrain and discriminate between 
different models of the physics of Type la SNe. We notice, in 
general, a larger dispersion in the U - B colour at all epochs 
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Fig. 6 ^-corrections to restframe Z?-band at redshift z = 0.3 (top 
panel) and z = 0.6 (bottom panel) using our templates (solid 
lines) and salt2 templates (dashed lines). 



than in all other colours. Due to the few data available, it is very 
difficult to establish whether the larger dispersion is a real char- 
acteristic of the £/-band or a consequence of systematic effects. 
The U -band photometry is affected generally by larger measure- 
ment uncertainties, smaller instrumental sensitivity, and larger 
extinction corrections (both atmospheric extinction, and inter- 
stellar extinction in the Milky Way and in the host galaxy). The 
broad variability of the effective f/-band bandpass shapes at dif- 
ferent telescopes could lead to larger dispersio n if accurate S- 
correc tions are not applied (see e.g appendix in Stan ishev et al.l 
(I2007ln . Even larger are perhaps the uncertainties that could be 
introduced by the colour term in the determination of the zero 
point. This is exp ected to be of the order of 0. 1 mag in the U- 
band photometry dSuntzeffll2000l) . Many of these systematic un- 
certainties are correlated at different epochs, and thus, can easily 
simulate a correlated intrinsic dispersion as the one observed in 



25 



20 



1 5 



10 



5 - 



0.0 



0.2 0.4 
E(B-V) 



0.6 



Fig. 7 Distribution of colour excess E(B - V) for the whole sam- 
ple including 80 SNe. 



i " " i i 1 1 1 1 

V This work Rv= 1 .00 
This work Rv=1.75 
CCM Rv = 3.1 
SALT2 




Fig. 8 The extinction law determined using the c param e ters in 
our model, compared with the one from Cardel li et al.l ( 1989) 
with Ry = 1.01 (light grey solid line) and Ry = 3.1 (dotted 
line). Also plot is the result by salt2 (dashed red line). The CCM 
law with Ry - 1.01 give the best fit to our data (solid points). 
The lighted shaded region represent the variability on the Ry 
parameter obtained by running the Monte Carlo simulation. The 
darker shaded region is obtained by adding 1 1 highly extincted 
SNe, and the CCM law with R v = 1 .75 give the best fit to our 
data (solid squares). A value of E(B - V) — 0.1 was assumed for 
all the curves in the plot. 



the data. This is certainly a strong concern for [/-band, but it 
could be a second order effect for other bands as well. A larger 
and more homogeneous data sample, such as the one collected 
by the SDSS collaboration, would be a major improvement for 
investigating these hypotheses. 

We have built colour curves for hsiao and salt2 spectral tem- 
plates by computing synthetic photometry, and compared them 
with our results. The left-hand panel of Fig.|9]shows a compari- 
son of the colour curves U-B, B-V,V-R and R-I for an s — 1 
supernova as derived from the three spectral templates consid- 
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ered. The right-hand panel shows the colour curves for differ- 
ent stretch values for salt2. We note that there is a general good 
agreement in B - V colour curves for a standard s = 1 supernova, 
with differences within the intrinsic colour dispersion. This is 
expected, since the availability of both B and V band data is ex- 
cellent. The larger deviations between models are present in the 
U - B colour curve. We note in particular the strong dependence 
on stretch shown by salt2 around day +40. This is unexpected 
and most probably an artifact of a poorly constrained model due 
too sparse data in the {/-band at this late epochs. The same rea- 
soning holds for the behaviour at early epochs. The salt2 colour 
curves show overall more wiggles than the colour curves derived 
in this work (see Fig. |2j». One possible cause is the inhomogene- 
ity of data quality. Some of the supernovae in the sample have 
very small measurement errors and if no intrinsic color disper- 
sion is considered, they weight very heavily in the overall fit, 
possibly biasing the model. The comparison between spectral 
templates shown in Fig. [5] leads to the same conclusions. The 
salt2 templates show a behaviour below 4000 A which is hard 
to attribute to real spectral features, both for early and late time 
(epoch=-10 and epoch=+20). These fluctuations are not present 
in the hsiao templates. Another anomaly in the salt2 templates is 
the emission feature at ~ 7600 A, shown in many of the spectra 
after maximum. Note, however, that these anomalies are within 
the reported uncertainties in the salt2 model. 

When comparing with other published results, one should 
keep in mind the different aims of the various analyses. The 
salt2 templates, for instance, are a side product of a light curve 
fitter. Thus, the colour curves derived by computing synthetic 
photometry are not built on colour data, but are a consequence 
of adjusting the spectra by the colours of the supernovae. In this 
sense, they are subject to uncertainties introduced by the tech- 
nique used for warping the spectra. Similarly the hsiao templates 
are built with the only aim of computing ^-corrections, and not 
for fitting lightcurves or exploring colour properties. As men- 
tioned in section [5] the averag e colour imposed to this template 
is the one by Knop et al. (2003) which based the U-B colour on 
the measurements of five supernovae at maximum. Another im- 
portant difference is the much larger number of parameters used 
for building salt2 spectral templates, leading to increased un- 
certainties compared to this work and the hsiao templates. When 
using few parameters, one assumes implicitly some "resonable" 
constrains on the models, such as smoothness and the number of 
knots for the spline interpolation. The ability to relax these as- 
sumptions is limited by the quality and ammount of the avalaible 
data. 

In section [6] we derive a reddening law from our sample 
which is well described by the CCM law with a low value of 
R v = 1.01 ± 0.25 for SNe with E(B - V) < 0.25 and R v = 1.75 
for SNe with E(B - V) < 0.7 mag. It should be noted that, we 
have not used the CCM law in any part of the analysis, except for 
correcting the data for galactic reddening. The spectral templates 
are warped to match the average colour using a spline interpola- 
tion, and not the CCM curve. We assumed no prior on any value 
of Ry or on any behaviour of the colour excess E(B - V). Yet, 
we find a good match to the CCM law. The direct comparison 
with what found by salt2, shows a very good agreement, except 
for the [/-band (see Fig. O, with differences of about 0.06 mag, 
for E(B —V) = 0.1. Given the intrinsic dispersion in the U-B, 
reported in Table [2]this is not a large discrepancy. Moreover, be- 
sides all the concerns already discussed about possible problems 
in the {/-band for the salt2 templates, one should also consid- 



ered the additional uncertainties specific to {/-band discussed 
above. 

When adding the 1 1 highly reddened SNe in our sample, up 
to E(B-V) < 0.7 ma^fl, the derived extinction properties change, 
yielding a larger value of the total-to-selective extinction param- 
eter Ry = 1.75 ± 0.27. The two values are statistically incom- 
patible, given that they are determined using largely overlapping 
samples. Moreover, as pointed out in section|7] the quoted uncer- 
tainty may be overestimated, since determined assuming uncor- 
rected intrinsic dispersion between SN epochs. A lower uncer- 
tainty is found for the fully correlated case, cr Rv ~ 0.13, which 
would make the difference even more significant. 

Using highly extincted SNe gives more leverage when study- 
ing extinction law properties. As it is not possible to perform 
this analysis only on the very small highly extincted SN sample 
available, the determination of Ry on the larger sample remains 
our most robust global estimate of the total to selective extinc- 
tion parameter. 

Based on our findings we argue that it may not be surpris- 
ing that different analyses, based on overall minimization tech- 
niques of the Hubble diagram, report different global fits of Ry, 
as its determination depends much on the selection of the sam- 
ple, and on the w ay the uncertaint i es are estimated. It is inter- 
esting to note that Tripp & Branch] d 19991) found a very similar 
result when minimizing the dispersion in the Hubble diagram: 
Ry = 0.876 (inferred from the measured Rb - 1.876 when con- 
sidering SNe with E(B - V) < 0.20 mag and Ry = 1.439 (in- 
ferred from the measured Rb - 2.439) when this condition was 
relaxed. 

As most samples used for fitting cosmological parameters, 
are selected based on low reddening, caution should be used 
when choosing the value of Ry for reddening corrections. Even 
though the corrections could be small, a systematic bias may be 
introduced, and could lead to substantially different results, e.g. 
as shown bv lConlev et al.l d2007l) . 

The many studies reporting measurements of Ry inconsistent 
with the Milky-Way value {Ry = 3.1) underline the need for a 
better understanding of dust properties in other galaxies. In this 
work we find that the reddening law for the low-extinction sam- 
ple is statistically incompatible with the value derived with the 
full available sample. This may be taken as suggestive evidence 
for additional processes involved in forming the effective red- 
dening law of SNe la besides dimming by interstellar dust, e.g. 
scatte red light echoes due to dust in the circumstella r environ- 
ment (IWand2005t iPatat et al.l2006l IWang et al.l20 07). or intrin- 
sic properties of the SN explosion mimicking reddening by dust. 
Unfortunatelly, only few high-quality colour measurements are 
currently available for highly reddened SNe la. Additional ob- 
servations of SNe la with E(B - V) > 0.2 should be carried out 
to further investigate this issue. 
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7 We note that a few SNe la that show even higher reddening have 
recently been observed ( Elias-Rosa et al . 2006, 2008). The values of 
R v determined for these individual SNe are in agreement with the larger 
value obtained in our work. 
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Fig. 9 Comparison of synthetic colour curves for all templates for s — 1 supernova (left panel). Synthetic colour curve evolution for 
different value of the stretch factor for the salt2 templates (right panel). 



Table 5 U — B: time evolution of the a and b parameters. cjj_ b = 
0.82 ± 0.06 



Epoch 


bit) 


a 


(t) 


-9 


-0.614 ± 


0.011 


-1.046 


± 0.101 


-6 


-0.594 ± 


0.011 


-1.002 


± 0.090 


-3 


-0.563 ± 


0.010 


-0.940 


± 0.077 





-0.462 ± 


0.008 


-0.779 


± 0.058 


3 


-0.325 + 


0.009 


-0.635 


± 0.078 


6 


-0.196 + 


0.009 


-0.648 


± 0.093 


9 


-0.102 ± 


0.008 


-0.848 


± 0.097 


12 


-0.032 ± 


0.009 


-1.074 


±0.125 


15 


0.025 ± 


0.011 


-1.153 


±0.145 


18 


0.079 ± 


0.011 


-0.959 


±0.136 


21 


0.126 ± 


0.011 


-0.586 


±0.151 


24 


0.161 ± 


0.012 


-0.183 


±0.184 


27 


0.173 ± 


0.012 


0.107 


±0.174 


30 


0.165 ± 


0.011 


0.260 ±0.136 


33 


0.149 ± 


0.013 


0.311 


±0.124 


36 


0.136 ± 


0.014 


0.305 


±0.134 


39 


0.128 ± 


0.014 


0.274 


± 0.128 


42 


0.123 ± 


0.012 


0.241 


±0.129 


45 


0.120 ± 


0.012 


0.224 


±0.135 


48 


0.118 ± 


0.014 


0.217 


±0.140 


51 


0.115 ± 


0.018 


0.209 


±0.145 


54 


0.113 ± 


0.023 


0.202 


±0.151 


57 


0.111 ± 


0.028 


0.194 


±0.158 


60 


0.109 ± 


0.033 


0.187 


±0.166 


63 


0.107 ± 


0.039 


0.180 


± 0.174 



Table 6 B - V : time evolution of the a and b parameters, cb-v - 1 
by construction. 



Epoch 


Kt) 


a 


(0 


-9 


-0.116 ± 


0.006 


-0.207 


± 0.050 


-6 


-0.104 ± 


0.005 


-0.187 


± 0.045 


-3 


-0.092 ± 


0.004 


-0.167 


± 0.040 





-0.067 ± 


0.004 


-0.135 


± 0.033 


3 


0.001 ± 


0.004 


-0.116 


± 0.034 


6 


0.089 ± 


0.005 


-0.267 


± 0.041 


9 


0.186 ± 


0.004 


-0.626 


± 0.040 


12 


0.306 ± 


0.005 


-1.033 


± 0.052 


15 


0.462 ± 


0.006 


-1.316 


± 0.058 


18 


0.652 ± 


0.005 


-1.380 


±0.051 


21 


0.838 ± 


0.005 


-1.272 


± 0.056 


24 


0.986 ± 


0.006 


-1.038 


± 0.064 


27 


1.067 ± 


0.006 


-0.727 


± 0.054 


30 


1.097 ± 


0.006 


-0.398 


± 0.046 


33 


1.092 ± 


0.008 


-0.109 


± 0.059 


36 


1.071 ± 


0.008 


0.082 


± 0.066 


39 


1.044 ± 


0.007 


0.156 


± 0.052 


42 


1.014 ± 


0.008 


0.158 


± 0.030 


45 


0.980 ± 


0.008 


0.134 


± 0.016 


48 


0.944 ± 


0.007 


0.122 


± 0.014 


51 


0.906 ± 


0.007 


0.123 


± 0.014 


54 


0.868 ± 


0.009 


0.128 


± 0.019 


57 


0.830 ± 


0.013 


0.129 


± 0.021 


60 


0.791 ± 


0.018 


0.131 


± 0.023 


63 


0.752 ± 


0.023 


0.133 


± 0.025 
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Table 7 V - R: time evolution of the a and b parameters. c V -r = 
0.37 + 0.03 





h(f\ 


a(f\ 


_q 

y 


ft 1 ft + ft ftf)7 

U.U1U in U.UU/ 


_n sat. + n 04.Q 
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o 


ft ftftS -+- ft ftftft 
U.UUO It u.uuo 


f) ^ 1 T, + ft D44 

U.J1J It U.UH-H 


j 


ft ftf)ft + ft rift's 
u.uuu in u.uuj 


-ft 4.84. + n rnQ 

U.tOt It U.UJ7 


n 


ft ftftft + ft ftft"S 


-ft 41 x + n rnft 

U.tlO It U.UJU 


J 


-ft f)9ft + ft ftf)7 
U.UiU it u.uu/ 


-ft 947 + ft 0S0 

\J.zJ-rl It U.UJU 


U 


-ft ITS8 + ft ftft7 
u.ujo m u.uu/ 


-ft Dft4 + ft OSd 

U.UUt It U.UJt 


q 

y 


-ft 1ft 1 + ft ftftft 


ft D1 8 + ft OSQ 
u.uio it \j.\jjy 


12 


-ft 1 1 8 + ft ftf)7 
U. 1 lO it u.uu/ 


_n 1 97 + ft ft77 
U.li/ it u.u/ / 


i j 


-ft f)89 + ft ftftft 
u.uoz, m u.uuo 


-ft ^ft^ + ft ftftl 

U.JUJ It U.UU1 


1 8 
l o 


ft ftft 1 + ft ftftS 
O.OUl It u.uuo 


fl qfn + n ft79 
U.7U/ It u.u/z. 


21 


ft Iftl + ft ftft7 
u.iui m u.uu/ 


-1 f)1ft + ft ftft4 
-1.U1U It u.uu^t 


94 


ft 1 Q7 + ft ftftS 

U. 17 / It U.UUO 


ft 898 -+- ft ftftS 

U.OZ.O It U.UUJ 


97 


ft 96? + ft DDR 
U.ZUi It u.uuo 


-ft ^74 + ft ftftft 
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Table 1. continued. 
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